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ABSTRACT 

We investigate the structure and evolution of a geometrically thin viscous Keplerian circumbinary 
(CB) disk, using detailed models of their radiative/convective vertical structure. We use a simplified 
description for the evolution of the cataclysmic binary and focus on cases where the circumbinary disk 
causes accelerated mass transfer (> 10~ 8 M 0 yr -1 ). The inner edge of the disk is assumed to be 
determined by the tidal truncation radius and the mass input rate into the disk is assumed to be a small 
fraction (10 _5 -0.01) of the mass transfer rate. Under the action of the viscous stresses in the disk the 
matter drifts outward with the optically thick region extending to several AU. The inner part of the 
disk is cool with maximum effective temperatures < 3,000 K while the outermost parts of the disk are 
< 30 K and optically thin. We calculate the effects of thermal instability on a sufficiently massive CB 
disk. It leads to outbursts reminiscent of those in thermally unstable accretion disks, with the instability 
remaining confined to the inner regions of the CB disk. However, for most of the evolutionary sequences 
the surface densities required to trigger instability are not reached. The spectral energy distributions 
from circumbinary disks are calculated, and the prospects for the detection of such disks in the infrared 
and submm wavelength regions are discussed. 

Subject headings: binaries: close — stars: evolution —- stars: cataclysmic variables 


1. INTRODUCTION 

Over the last decade evidence has been accumulating 
suggesting the presence of gaseous matter beyond the orbit 
of the components in cataclysmic variable binary systems 
(CVs). Its existence has been inferred from P-Cygni pro¬ 
files in nova-like variables, super soft sources, and dwarf 
novae in outburst (see, for example, Deufel et al. 1999). 
In addition, the observation of narrow widths of single- 
peaked emission lines in SW Sex stars (see Thorstensen et 
al. 1991; Hellier 2000) and low velocity spectral features 
in AM CVn which do not follow the orbital motion of the 
binary system (see Solheim & Sion 1994) provide further 
support for this interpretation. Such a gaseous compo¬ 
nent may take the form of an extended shell or a flattened 
distribution of matter surrounding the binary system. As¬ 
suming that some fraction of this matter remains bound to 
the system and condenses into the orbital plane, in anal¬ 
ogy to the compression of a radiation-driven wind into the 
equatorial plane in Be-stars (Bjorkman & Cassinelli 1993, 
Bjorkman & Wood 1995), a circumbinary disk may form. 
The presence of such a CB disk, as proposed by Spruit 
& Taam (2001), can provide an attractive explanation for 
some of the discrepancies between observations and the 
current theory of evolution of CVs. 

The importance of a CB disk for the evolution of a CV 
lies in the angular momentum it can extract from the bi¬ 
nary, leading to mass transfer rates exceeding those due 
to magnetic braking and/or gravitational radiation. The 
binary exerts gravitational torques on the external disk, 
thereby transferring some of its orbital angular momen¬ 
tum. If some fraction of the mass lost from the main 
sequence-like secondary of the system forms a CB disk, a 


feedback process can operate to accelerate the binary evo¬ 
lution. As shown in the detailed computations by Taam 
& Spruit (2001) the evolution of the system can be sig¬ 
nificantly accelerated and the mass transfer rate from the 
secondary towards the white dwarf, M 2 , elevated to rates 
as high as 10~' M 0 yr -1 . 

In previous studies, we have focused on the response of 
the mass losing component and on the evolution of the bi¬ 
nary system resulting from the influence of the CB disk, 
using a one zone approximation for its vertical structure. 
Here, we report on the results of detailed numerical cal¬ 
culations for the vertical structure and radial evolution of 
such disks surrounding a model CV. In order to determine 
the observability of such disks, their spectral energy distri¬ 
butions are also presented. The assumptions of the model 
computations are described in the next section. In §3 we 
describe the evolution of a CB disk neglecting angular mo¬ 
mentum losses from the binary system. Results with the 
binary evolution taken into account are shown in §4 and 
several spectral energy distributions for some illustrative 
systems are presented in §5. We summarize and conclude 
in the last section. 

2. MODEL 

We assume that the CB disk is axisymmetric and ge¬ 
ometrically thin such that matter rotates at the Ke¬ 
plerian speed. The formulation of the equations for 
mass, angular momentum, and energy transport is as de¬ 
scribed in Hameury et al. (1998). Models for the ra¬ 
diative/convective vertical structure have been calculated 
and tabulated in a grid for a range of radii, mid-plane 
temperatures, and column densities. The global structure 
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and evolution of the disk follows from the coupling of this 
grid to the time dependent transport of mass, energy, and 
angular momentum in the radial direction. The viscous 
stress tensor, r r ^, is assumed to be equal to aP where a is 
the viscosity parameter introduced by Shakura & Sunyaev 
(1973) and P is the local pressure in the disk. The calcula¬ 
tions of the disk differ from those of Taam & Spruit (2001) 
where the vertical structure of the disk was calculated in 
a one zone approximation. Furthermore, we take into ac¬ 
count departures from local thermal equilibrium and in¬ 
clude convective energy transport. The latter is treated in 
the mixing length approximation with the ratio of mixing 
length to pressure scale height taken to be equal to 1.5, 
where the functional form for the pressure scale height is 
chosen to ensure that it does not exceed the vertical scale 
height. 

Optically thin regions above the disk photosphere are 
included in the computation and are treated in the grey 
approximation. The input physics includes an equation of 
state of matter in which the ionization/dissociation bal¬ 
ance of hydrogen is taken into account. The opacities are 
taken from Bell et al. (1997) and include contributions 
from iron poor opacities at temperatures (T < 1500 K) 
where dust particles exist (Henning & Stognienko 1996). 
It is assumed that the composition of the material corre¬ 
sponds to a solar mix ( X = 0.7, Y = 0.28, and Z = 0.02). 

The boundary conditions are formulated such that both 
the inner radius and outer radius of the disk can vary in 
time. At the inner radius, it is assumed that the mass 
input rate into the CB disk is a constant fraction of the 
mass transfer rate from the secondary M\ = — 8M 2 . Al¬ 
though the mass input into the disk may extend over a 
range of radii, we have adopted the simplest approach and 
assumed that the mass is deposited at its innermost edge. 
The condition that this input mass flux is carried outward 
by viscous diffusion serves as an inner boundary condition, 
which is formulated as a relation between S and its gradi¬ 
ent. The specific angular momentum of the mass fed into 
the disk is taken to be equal to that of a Kepler orbit at 
the inner edge. The mass deposited is thus neutral in an 
angular momentum sense. In addition, it is assumed that 
the radial temperature gradient vanishes. Our treatment 
differs from that adopted in Pringle (1991) in that we have 
a mass input into the disk which is absent in the formula¬ 
tion developed by Pringle (1991). In the latter equivalent 
approach, the radial velocity at the inner edge is assumed 
to vanish. The column density and the mass outflow rate 
at the outer boundary are set to 10~ 2 g cm -2 and zero re¬ 
spectively. With these treatments of the boundaries, the 
disk structure and evolutionary equations are solved in a 
fully implicit fashion on an adaptive grid (typically with 
400 zones). For a full description of the numerical method, 
see Hameury et al. (1998). 

Since the primary focus of this investigation is on the 
structure of the CB disk, the evolution of the binary sys¬ 
tem is treated in a simplified approximation wherein the 
response of the mass losing star is modeled assuming ho¬ 
mology (Spruit & Ritter 1983; see also Spruit & Taam 
2001). The mass transfer process is assumed to be con¬ 
servative. That is, a very small fraction, 6 < 0.01, of the 
mass lost from the donor star, which itself is assumed to be 
small compared to the mass transferred to its companion, 


is assumed to be deposited into the CB disk. 

Angular momentum loss associated with magnetic brak¬ 
ing is assumed to take place on a fixed timescale f w . The 
contribution of the CB disk to the angular momentum loss 
from the binary is given by the viscous torque at the inner 
edge of the disk, 

j d = 37r(i?i/a) 1/2 noa 2 ^iSi, (1) 

where flo is the orbital frequency of the binary system, a is 
the orbital separation of the binary components, and Rj, 
v.i, and are the radius, viscosity, and column density at 
the inner edge of the CB disk respectively. The evolution 
of the system follows from the total angular momentum 
loss, (see Spruit & Taam 2001): 

- J/J = t - 1 + 3 7 r(i? i /a) 1 / 2 (l + q)/M 2 v^ (2) 

It is assumed that the gravitational torque is concentrated 
near the inner edge of the disk and that the radius of the 
inner edge is approximately equal to 1.7 times the orbital 
separation. This latter assumption is based on the work 
of Artymowicz & Lubow (1994) who show that the ratio 
of the inner radius of the CB disk to the orbital separation 
is relatively insensitive to the mass ratio of the binary sys¬ 
tem, q, due to the steep decrease of the tidal torques with 
distance. 

3. CB DISKS IN NON-EVOLVING SYSTEMS 

As a first step we computed several CB disk models in 
which we neglected all the angular momentum losses from 
the binary system. The binary evolution is thereby frozen 
with the inner radius Ri and mass input rate M, (instead 
of 5) set as constant parameters. Our calculations start 
with a very low density (£ = 0.01 g cm -2 ) ring of matter 
at the inner edge, and we follow the expansion of the disk 
as matter gradually piles up and diffuses outwards. 

Fig. 1 is an example of the evolution of the column den¬ 
sity at the inner radius, Si, and of the outer radius, Rd, for 
a CB disk with a — 0.001, R-, = 10 11 cm, M\ = —10 14 g s -1 
(the negative value indicating mass is flowing away from 
the binary) and a binary mass M of 1 M 0 . The outer 
disk radius as shown in Fig. 1 is defined as the radius 
at which the disk becomes optically thin (r = 1). After a 
transitory period lasting at most a few hundred years (vis¬ 
cous timescale of the inner edge), both Si and Rd increase 
gradually in a roughly self-similar, power law-like fashion. 
This behaviour is independent of the initial setup and, to 
a large extent, independent of the outer boundary condi¬ 
tions on M and S. Increasing the mass input rate speeds 
up the evolution with higher values of S; reached sooner. 
Changing a has a weaker influence with lower viscosities 
leading to slightly higher values of Si. Both of these effects 
can be seen in Figs. 3-6 which take the binary evolution 
into account and are discussed in the next section. 

The evolution of the column density, mid-plane temper¬ 
ature, effective temperature, ratio of scale height to radius, 
and mass flow rate are illustrated as a function of radius 
and time in Fig. 2. It is evident that the column densities 
continually increase at a given radius and reach densities 
as high as 10 4 g cm~ 2 at the inner disk edge. The point at 
which the optical depth equals unity (denoted as a solid 
dot), corresponding to Rd in Fig. 1, also increases with 
time and the optically thick disk reaches a size of about 10 
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AU at the end of the calculation. The mid-plane temper¬ 
atures (solid line) and the effective temperatures (dashed 
line) both increase with time, reflecting the increasing sur¬ 
face density and energy dissipation in the disk associated 
with the increasing mass flow rate. The mid-plane (effec¬ 
tive) temperatures decrease monotonically with distance, 
from 10000 (3000) K in the innermost regions to 20 (10) K 
where the disk becomes optically thin. However, we note 
that such densities and temperatures are reached only af¬ 
ter a Hubble time, at which point the neglect of angular 
momentum losses from the binary is not realistic. The 
disk remains geometrically thin throughout the evolution, 
with the ratio of pressure scale height to radius in the op¬ 
tically thick region, H/R, < 0.05. We find that H/R does 
not increase uniformly with radius, so that the outer disk 
region is in the shadow of the inner region, reminiscent of 
the protostellar circumstellar disks studied by Bell et al. 
(1997). 

The simple self-similar CB disk model of Spruit & Taam 
(2001), which assumes that the viscosity is a linear func¬ 
tion of R only, v = lrfR/R/), gives for the inner regions of 
the disk: 


E(R,t) 


t \ 1/2 (-M;) / R \ 3/2 
tyi ) ZlTVi \Rij 


( 3 ) 


where t v ; = R? /is the viscous timescale at the inner edge 
of the disk. This implies a scaling for E, oc M\/ ffa, which 
we find to be roughly reproduced by the numerical results 
(see e.g. the early evolution of the disks in Figs. 5-6). 


4. CB DISKS IN EVOLVING SYSTEMS 

Evolutionary sequences were calculated for model CV 
systems consisting of a secondary component of 0.55 M 0 
and a primary white dwarf component of 0.95 M 0 . The 
disk calculations were performed for several fractional 
mass input rates into the CB disk 8 < 0.01, viscosity pa¬ 
rameters a < 0.01, and time scales for angular momentum 
loss by magnetic braking t w < 10 10 yrs. 

4.1. Acceleration of the binary evolution 

The results of an evolutionary sequence with parameter 
values 6 = 0.005, t w = 2 x 10 8 years, and a = 0.001 are 
illustrated in Figs. 3 and 4. The CB disk initially evolves 
in the same fashion as described in the previous section. 
Angular momentum losses by magnetic braking then start 
changing the system parameters and the system evolves 
in a standard way, with —M; decreasing together with the 
mass of the secondary. The CB torque begins to domi¬ 
nate the evolution at t ss 20 Myr, when the mass input 
rate into the CB disk is approximately at its minimum, 
Mi « -3 x 10" 11 M 0 yr- 1 (M 2 « 6 x 10" 9 M 0 yr" 1 ). 
The mass input rate is now proportional to the surface 
density since M-, /M 2 ~ J-JJ ~ — (i/E)i/M 2 , leading to a 
runaway on a timescale comparable to that on which the 
CB torque becomes dominant (Spruit & Taam 2001). The 
acceleration is manifest in Fig. 3 as an abrupt increase of 
—Mi to > 10 -9 M 0 yr -1 , leading to an increase in the or¬ 
bital period as a result of the tendency for the secondary 
to expand as it departs significantly from thermal equilib¬ 
rium. 

By the end of the calculation at t « 65 Myr the orbital 
period is ~ 20 hr and the secondary has M 2 « 0.03 M 0 . 


The response of the star to mass loss certainly deviates 
from homology during the runaway but there is little doubt 
that it will end with the dissolution of the secondary. Qual¬ 
itatively similar results were obtained by Taam & Spruit 
( 2001 ) . 


4.2. Influence of the different parameters 

The influence of changing a is also shown in Fig. 3. The 
CB torque becomes dominant when Jd ~ J w - Combining 
Eq. 1-3, we find this happens when (Spruit & Taam 2001) 

Jd/J ~ (tAvi) 1/2 Mi/M 2 - -f" 1 (4) 

For a given magnetic braking torque and Mi, the disk 
torque dominates at t oc f v ; oc 1/a. This scaling is con¬ 
sistent with the numerical results of Fig. 3. Note that for 
very low values of a (< 0.0001 for <5 = 0.005 and t w = 0.2 
Gyr) the CB torque would be entirely negligible during 
the course of the evolution. It should be pointed out that 
we have implicitly assumed that a is constant throughout 
the disk. It is possible that a may decrease significantly at 
large radii (see Gammie & Menou 1998) where the mag¬ 
netic Reynolds number is decreased below values neces¬ 
sary for a magnetic dynamo to operate. In this case, the 
size and the properties of the outer disk would be modi¬ 
fied from that calculated here. However, the extent of the 
magnetically decoupled region is model dependent as it is 
very sensitive to the ionization fraction in the disk (see 
Fromang, Terquem, & Balbus 2002). 

Changing the magnetic braking timescale t w changes M\ 
since the initial rate is set by M 2 /M 2 ~ J/J rss tfff (Eq. 2). 
Using this in Eq. 4, we find the CB torque should thus start 
dominating the evolution at roughly the same time regard¬ 
less of t w . This is confirmed by the calculations shown in 
Fig. 5 in which the disk torque dominates after t « 20 Myrs 
for f w = 0.2, 1 and 5 Gyrs. Note that E; at t « 20 Myrs 
decreases with increasing f w since E, oc —Mi oc f ” 1 during 
the early magnetic braking-dominated evolution. The con¬ 
ditions for the runaway should be independent of t w and 
we indeed find the disks in Fig. 5 to have similar properties 
when it occurs (E, w 40000 g cm -2 , R « 10 AU). 

Fig. 6 illustrates the influence of changing the fraction 
of the mass transfered to the CB disk. For (5 > 0.01 the 
system dissolves in < 1 Myr suggesting that systems with 
high values of 6 are unlikely to be observed. The timescale 
on which the disk torque dominates dramatically increases 
as 8 becomes smaller. The effect is much stronger than de¬ 
creasing a as already noted by Taam & Spruit (2001). For 
values of 8 < 0.0001 the disk torque becomes important on 
timescales longer or comparable to the magnetic braking 
timescale f w . The secondary has then already lost most of 
its mass so that the CB disk has no impact on the binary 
evolution. Being smaller and less massive, the CB disk is 
also fainter and more difficult to detect (see §5). 

Taam & Spruit (2001), using a detailed evolutionary 
code for the secondary but a simpler disk model, found 
that the CB torque was negligible when 8 < 0.01. The 
present calculations indicate that disks formed with 8 > 
0.0001 can still influence the binary. The discrepancy can 
be traced to the assumptions made in the different models. 
In Taam & Spruit (2001), the column densities at the inner 
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edge of the disk are very close to the lower turning point 
of the S curve (E max ~ 4000 g cm -2 at R = 10 11 cm with 
M — 1 M 0 and a = 0.001, see Fig. 1 of Taarn & Spruit). 
The specific interpolation scheme between the two stable 
branches which they used forces the column density to 
stagnate around E max once this value is reached. The de¬ 
tailed vertical structure computations, which include con¬ 
vection, show that the actual value of E max is an order 
of magnitude higher than in the one zone radiative model 
(E max « 40000 g cm -2 with the same parameters as above 
using Eq. 32 of Hameury et al. 1998 ; see Fig. 2 of Lasota 
2001 for a discussion of the importance of convection on 
Emax)- Much higher densities can be reached on the cold 
branch and hence higher torques can be achieved during 
the evolution for the same <5 and a. On the other hand, 
the approximate treatment for the response of the donor 
to mass loss for the binary evolution in the present calcu¬ 
lation may artificially decrease the minimum S needed to 
have a dominant CB torque. 

4.3. Can CB disks become thermally unstable ? 

Depending upon the evolution of the CB disk, E; may 
become higher than the critical column density, E max , 
above which a cold, optically thick and geometrically thin 
disk becomes thermally and viscously unstable. The insta¬ 
bility is due to the large change in opacity which accom¬ 
panies hydrogen ionization (see Lasota 2001 for a review). 
Such a CB disk would present limit cycle oscillations anal¬ 
ogous to those which occur in the accretion disks of dwarf 
novae. 

While the linear stability problem is the same for the 
accretion disk and the CB disk, one should expect some 
differences to show up in the nonlinear development. In 
the accreting case, the increased mass flux in the hot state 
propagates inward into a decreasing volume, amplifying 
its effects. The instability in a CB disk starts at the inner 
edge and propagates outward into an increasing volume. 
One might thus expect its effects to be milder, and the 
transition front to stall at some moderate distance from 
the inner edge. To check on these expectations we have 
followed the instability in a few cases. But first, we inves¬ 
tigate for which parameters the evolving binary with a CB 
disk can reach instability. 

We combine the expression for E(i?;, t) (Eq. 3) with that 
for E max (Eq. 32 of Hameury et al. 1998) to estimate the 
time needed to reach this critical density, assuming the 
binary evolution is frozen (constant Ri and Mi). The cen¬ 
tral temperature at the critical density is T c ss 20000 K 
(Hameury et al. 1998). We use this temperature to get an 
estimate of « q:c 2 /Ok, where c 2 and Ok are respectively 
the sound speed and Keplerian frequency at the inner edge. 
We find 

t~ 6 x 10 10 Mf 12 M;" 2 Rf^ years (5) 

where a is in units of 0.001, M; is in units of —10 14 g s _1 , 
M is in solar masses and R\ is in units of 10 11 cm. The 
results from a few numerical tests show a shallower depen¬ 
dence on M[ (-1 instead of -2) but this is still a fair order- 
of-magnitude approximation. The steep dependence on 
Ri implies that instability can be reached within a Hubble 
time only if the binary is sufficiently tight during most of 
its evolution. 


In the calculations of Figs. 3-6, the CB disk never 
reaches the critical densities required to become unsta¬ 
ble. For reasonable assumptions on M\ (0.6-1.0 M 0 ) 
and M 2 (0.1-0.6 M 0 ) the CB disk has an initial Ri of 
ss 0.01 AU for which the critical density E max is of the 
order of 10 5 g cm' 2 for a = 0.001. The density in the 
disk does not build up fast enough to trigger the instabil¬ 
ity before the binary starts evolving away from the initial 
conditions. Unrealistically high values of S > 0.1 would be 
required to quicken the density buildup. At these rates, 
the binary would also dissolve in a very short time. 

If magnetic braking dominates during the evolution, the 
CB disk is necessarily not very massive and E; stays far 
from its critical value. If the disk torque dominates, the 
mass transfer rate and Ei increase quickly, particularly 
during the runaway. However, under the present assump¬ 
tions, the secondary expands as it rapidly sheds mass onto 
the white dwarf (Eq. 14 of Spruit & Ritter 1983). The 
runaway is thus accompanied by an increase of the orbital 
period, and hence of R\ (see Fig. 3), which prevents Ei 
from becoming greater than E max (oc Ri). 

The adiabatic mass-radius exponent is -1/3 for a se¬ 
quence of homologous stars; in this respect they behave 
like convective stars (Spruit & Ritter 1983). This limits 
the freedom of the model. We have relaxed this in a test 
calculation in which the mass-radius coefficient of -1/3 as it 
appears in the model was increased to +1/2. The other pa¬ 
rameters are a = 0.005 and 5 = 0.005 (Fig. 7). The global 
evolution is very similar to the one shown in Fig. 4 except 
that the secondary contracts in response to the runaway 
mass transfer at t « 30 Myr. The orbital period decreases 
to ~ 2 hours, bringing the inner edge of the CB disk to 
0.007 AU (M 2 « 0.085 M 0 , R 2 « 0.16 R 0 ). The inner 
column density E; is then greater than the critical den¬ 
sity and the CB disk, which is now ~ 0.003 M 0 , becomes 
unstable. 

We assume that, when the CB disk is hot, the viscosity 
coefficient a^ot increases to 0.01 by analogy with the mod¬ 
els of dwarf novae outbursts. The inner region of the CB 
disk cycles between a hot and a cold state in a complex 
fashion on a timescale of ~ 15 years, with the hot state 
lasting about a year (see Fig. 7). Changes in the CB disk 
on such a short time scale do not translate into the mass 
transfer rate, which changes only on time scales of the or¬ 
der the ratio of the scale height of the stellar atmosphere 
to its radius times t w . The mass transfer rate is therefore 
kept fixed at its current value when the instability starts 
(M 2 « 2.4 x 10~ 7 M 0 yr _1 ). 

The instability would therefore only show up, if at all, 
as a brightening of the system in the near-infrared. The 
CB disk can indeed become very bright with a bolomet- 
ric luminosity increasing by almost an order of magnitude 
to about 2 x 10 34 erg s _1 . The bolometric luminosity is 
calculated assuming the optically thick regions of the disk 
radiate as a blackbody and that the disk is seen face-on. 
Fig. 8 shows the evolution of the CB disk during the rise 
to the high luminosity state. A front develops at the in¬ 
ner unstable radius and propagates outward, transporting 
matter to larger radii and gradually heating the inner re¬ 
gions to mid-plane temperatures reaching 10 s K. This is 
analoguous to an outside-in outburst in a dwarf novae ac¬ 
cretion disk, where the outer (tidally truncated) radius 
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becomes unstable first and a front propagates in the direc¬ 
tion of the flow. However, in the CB case, the front stalls 
quickly due to the combination of the higher critical den¬ 
sity needed to heat the disk and lower available density to 
do this (dilution into a larger volume of the transported 
mass). In contrast, a dwarf novae outside-in front heats 
the entire accretion disk because of the decreasing E max 
and volume at smaller radii. At the luminosity maximum, 
only the inner 0.02 AU are hot compared to the total (op¬ 
tically thick) CB disk size of 2.5 AU. 

Our choice of boundary condition precludes the possi¬ 
bility of mass flow to the secondary. This may not be 
exact during outburst when the inner region becomes suf¬ 
ficiently thick ( H/R > 0.05) that some matter may flow 
from the CB disk back to the companion star (Artymow- 
icz & Lubow 1996). Such flows are non axisymmetric and 
their inclusion is beyond the scope of our simplified one 
dimensional radial description of the vertically integrated 
structure. Note that these flows are unlikely to change our 
conclusions on the possibility of having outbursts in CB 
disks. 

5. SPECTRAL ENERGY DISTRIBUTIONS 

Since the temperatures, densities, and size of the CB 
disk vary with time, the spectral appearance will depend 
upon the evolutionary stage of the binary system. The 
energy distributions for a CB disk viewed normal to its 
surface for three model sequences are illustrated in Figs. 
9-11 assuming a distance of 100 pc. The distributions are 
shown for the CB disk early in the evolution, prior to, and 
during the accelerated mass transfer phase. In these calcu¬ 
lations the disk surface is assumed to emit black body ra¬ 
diation at its effective temperature. This should be a rea¬ 
sonable approximation since the frequency dependence of 
the opacities (see Bell et al. 1997) show the CB disk to be 
optically thick to radiation at wavelengths < 350pm. The 
figures show the integrated flux over the optically thick 
parts of the disk. 

It is seen that the continuum peaks in the near-infrared 
at A ~ 3pm for the models sequences during the entire 
binary evolution. There is a tendency for the flux to in¬ 
crease for the more advanced evolutionary phases, when 
CB disks are larger and more massive. However, the flux 
does not increase at all wavelengths. The flux at A < 3pm 
can decrease even though the mass input rate into the CB 
disk increases (see Fig. 9). This is due to the increase of 
the binary orbital separation in the accelerated evolution 
phase, causing the radiating surface area to increase and 
the temperatures in the inner regions to decrease. These 
results taken together indicate that the spectral energy 
distributions are not sensitive to a and d in the range of 
0.001 - 0.005. 

To determine the observability of such disks, the flux 
contributions from the donor star and a steady accretion 
disk surrounding the white dwarf star must be taken into 
account. As an illustration we take an evolution with 
a = 0.005, S = 0.005, and t w = 2 x 10 8 yr at t = 5 x 10 6 
yr. At that evolutionary phase, a secondary of mass equal 
to 0.36 M 0 transfers matter to the white dwarf of mass 
equal to 1.14 M 0 at the rate of 1.3x 10 18 g s _1 . We assume 
that the outer radius of the accretion disk is 0.67 times the 
Roche lobe radius of the white dwarf. The luminosity and 
effective temperature of the secondary are 1.1 x 10 -3 L 0 


and 3500 K respectively, the values for a main sequence 
star of the same mass. Assuming a main sequence sec¬ 
ondary is likely to overestimate the contribution from the 
donor star, which is not in thermal equilibrium because of 
the mass transfer. In Fig. 12 the spectral energy distribu¬ 
tion from a steady state accretion disk and the low mass 
secondary component are included with the contribution 
from the CB disk. It is evident that the contribution of 
the CB disk increases with wavelength, and exceeds that 
due to the accretion disk and donor by a factor of about 
20 and 400 at 10 pm and 350pm respectively. For wave¬ 
lengths where the CB disk dominates the wavelength de¬ 
pendence is, F\ ~ A -2 ' 6 . The radius enclosing 90% of the 
total flux emitted by the CB disk is ~ 5 x 10 11 — 1.6 x 10 12 
cm at 10pm and 10 12 — 10 13 cm at 100pm. Hence, the 
cool outer regions of the CB disk at R > 5 x 10 11 cm 
provide the infrared and subnnn waveband signatures of 
its presence. We point out that the neglect of irradiation 
effects by the secondary is not expected to significantly 
modify this result since its importance in the outer parts 
of the disk (as inferred from the variation of H/R with re¬ 
spect to radius) is diminished by the shadowing effect by 
the innermost regions. We note, however, that the actual 
height of the disk differs from that estimated by the pres¬ 
sure scale height and irradiation effects of the outermost 
regions could be important (see Chiang & Goldreich 1997). 
In this case the infrared and sub mm fluxes would be in¬ 
creased and, hence, the fluxes estimated from the CB disk 
in the absence of irradiation effects should be considered 
a lower bound. 

6. CONCLUSIONS 

The evolution of a CB disk has been investigated for 
cases where it significantly accelerates the binary evolu¬ 
tion, such that the mass transfer rate onto the white dwarf 
can reach M > 10 -8 M 0 yr -1 , representative of nova-like 
variables and super soft sources. The column densities 
and temperatures in these disks are found to be similar to 
those in circumstellar (CS) disks surrounding young stel¬ 
lar objects. The column densities decrease with radius 
and reach maxima in the innermost regions of the disk 
with T, max ranging up to 4 x 10 4 g cm -2 . The midplane 
temperatures range from ~ 3500 K at the inner edge to 
less than about 30 K in the optically thin outer regions. 
In contrast to the CS disks, where accretion can take place 
steadily, the outflowing mass in CB disks is distinctly non¬ 
steady. The outer boundary of the CB disk, Rd, increases 
with time (roughly as Rd ~ t 1 / 2 ). We find that CB disks 
around CVs should extend to several AU, about 100 times 
smaller than CS disks. 

The numerical results show that the fractional mass in¬ 
put rate into the CB disk required to accelerate the evo¬ 
lution is lower than found in Taam & Spruit (2001). This 
is due to the effect of convection on the disk structure (ig¬ 
nored in Taam & Spruit), which causes the lower stable 
branch of the v — E relation to extend to much higher 
column densities. 

We find that the CB disks can become thermally un¬ 
stable, but only under rather special conditions, namely 
at high mass transfer rates combined with short orbital 
periods. Such conditions could be met in the short pe¬ 
riod double degenerate AM CVn systems or, in a different 
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context, in disks around Be stars (Lee et al. 1991). For 
normal CVs, we find that the CB disk is unlikely to become 
unstable because of the long timescale needed to reach the 
critical column density. However, depending upon the re¬ 
sponse of the secondary to mass loss, the CB disk could 
become unstable during the accelerated evolution. If the 
CB disk is unstable, a heat front develops and propagates 
outward from the inner radius. In contrast to the accre¬ 
tion disks in dwarf novae, the heating front is quenched 
quickly because of the difficulty in raising the lower col¬ 
umn densities encountered at larger radii to the increasing 
critical column densities needed to sustain propagation. 
The inner region cycles between the hot and cold states 
on a timescale of tens of years, inducing changes in the 
bolometric luminosity of a factor ~ 10. 

We find that the spectral energy distributions expected 
of CB disks in CVs can dominate the emission from the 
donor star and the accretion disk of the white dwarf at 
wavelengths A > 3/irn. At longer wavelengths the relative 
contribution from the CB disk to the total emission from 
the system increases. 

Typical flux levels for a nova-like variable for an assumed 
distance of 100 pc would be about 0.3 Jy at 10 /mr and 
a factor of 10 lower at 450/im. Although the distances to 
CV’s are not well known, the model calculations suggest 
that even in nova-like variables with distances order of 1 
kpc, detection of CBs in the 10 - 20 /cm region should be 
possible using the infrared detectors on the Gemini Ob¬ 
servatory (Telesco et al. 2001) and the Keck Observatory 
(Jones et al. 1998). Although the emission from these 
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CB disks becomes optically thin at 350 /tm, it is possi¬ 
ble that the longer wavelength radiation at 450/im may be 
detectable using the SCUBA instrument. 

In our interpretation the nova-like variables would be 
CVs with mass transfer accelerated by a CB disk, which 
makes them natural targets for searches for CB disks. 
However, the predicted emission from the CB disk before 
the accelerated phase is similar, for a substantial period of 
time, to that during the accelerated phase. The possibility 
for their detection is therefore not limited to the systems 
with high mass transfer rates. We speculate that the ‘extra 
infrared components’ discovered by Harrison et al. (2000) 
in several well studied dwarf novae may in fact be due to 
CB disks. These observations show that searches of CB 
disks surrounding bright cataclysmic variable systems are 
technically possible with current instrumentation. 
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years 

Fig. 1.— The evolution of a model CB disk with constant mass input rate ( M[ = — 10 14 g s —1 ) and inner radius (R[ = 10 11 cm), for 
viscosity parameter a = 0.001. Top panel: column density at the inner disk edge Si. Bottom panel: outer radius of the disk defined as the 
radius where the optical depth r = 1. The radial structure at the times indicated by the dots on these curves are presented in Fig. 2. The 
surface density required for disk-instability (S ma x ~6x 10 4 g cm -2 in this case) is not reached. 
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Fig. 2.— The radial structure of the CB disk presented in Fig. 1 at different times during its evolution. The column density, E, mid-plane 
temperature T c (solid line), effective temperature T e ff (dashed line), ratio of pressure scale height to radius, H/R , and mass flow rate M (in 
units of 10“ 12 Mq yr _1 , negative values indicate mass is flowing away from the binary) are illustrated as a function of radius and time (each 
radial curve corresponds to a dot in Fig. 1). The radius at which the optical depth equals 1 is marked by a dot on each curve. 
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Fig. 3. — Evolution of the CB disk with the evolution of and angular momentum loss from the binary system taken into account. Three 
calculations are shown in each panel for values of the viscosity parameter a = 0.01, 0.005 and 0.001. The calculations start with a secondary 
mass of 0.55 M© and a white dwarf of 0.95 M©. Viscosity parameter is a = 0.001, time scale for angular momentum loss by magnetic 
braking of the binary is £ w = 2 x 10 8 years, and mass fed into the CB disk —Mi is a fraction S = 0.005 of the mass transfer rate. From top 
to bottom: the mass input rate into the CB disk —Mi; the mass of the secondary and the mass of the disk (dashed line); the orbital period; 
the outer radius of the CB disk (at which the optical depth equals 1) and the inner radius of the CB disk (dashed line); the column density 
at the inner disk edge Ei and the corresponding critical column density E ma x above which local thermal instabilities are expected (dashed 
line). The radial structure of the CB disk is shown at different times (indicated by the dots in the first panel) in Fig. 4. 
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Fig. 4.— The radial structure of the CB disk with a = 0.001 presented in Fig. 3, at different times during its evolution. The column density 
£, the mid-plane temperature T c (solid line) and the effective temperature T e ff (dashed line), the ratio of the pressure scale height to radius 
H/R , and the mass flow rate M (in units of 10 -10 Mq yr _1 ) are shown as functions of radius and time (at t = 5 x [10 5 ,10 6 ,10 7 ] years, see 
Fig. 3). The radius at which the optical depth r = 1 is shown as a dot on each curve. 
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Fig. 5.— Same as Fig. 3 but showing the influence of varying t w , the timescale for binary angular momentum loss by magnetic braking. 
The other parameters are a = 0.001 and 5 = 0.005 and the calculations start with a secondary mass of 0.55 M© and a white dwarf of 0.95 
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Fig. 6.— Same as Fig. 3 but showing the influence of varying <5, the fraction of the binary mass transfer rate which is input in the CB disk. 
The other parameters are a = 0.005 and i w = 2x 10® years and the calculations start with a secondary mass of 0.55 Mg and a white dwarf 
of 0.95 M m . 
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Fig. 7.— Example of a thermally and viscously unstable CB disk. The calculation starts with a 0.55 M© secondary and a 0.95 M© primary. 
We assume the secondary contracts in response to mass transfer (see §4.3), Cfccold = 0.005, o;hot = 0-01 and S = 0.005. The evolution up to 
t ~ 30 Myr is comparable to the model shown in Fig. 3. During the runaway, the radius of the secondary shrinks to 0.16 R©, decreasing 
the orbital period to ~ 2 hours and bringing the CB disk to shorter radii (0.007 AU). The critical E ma x is then low enough for the disk 
to become unstable. The upper and lower panels show the bolometric luminosity and the column density at the inner edge of the CB disk 
after it becomes unstable. The mass input into the disk is constant during the cycles at —M[ ^ 10 —9 M© yr —1 . The evolution of the radial 
structure of the CB disk for the outburst at t & 24 yrs is shown in Fig. 8. 
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Fig. 8. — Evolution of the radial structure of a CB disk during the rise to a luminous state. The model shown is that of Fig. 7 for the 
outburst peaking at t « 24 yrs. The top panel is a zoom on the bolometric luminosity lightcurve with the dots indicating at which times 
the 6 radial structures shown underneath were taken. The bottom three panels show the radial profiles of the column density E, mid-plane 
temperature T c and mass flow rate in the CB disk M (in units of 10 -6 Mq yr —1 ). The two dashed lines in column density panel are the 
critical densities E ma x (top) and E m i n (bottom). In contrast to the outbursts in standard accretion disks, the front does not propagate very 
far out into the CB disk before it stalls. 





15 



log X (micron) 


Fig. 9.— The spectral energy distribution for a CB disk for the model sequence with a = 0.001 and 5 = 0.005 for an assumed distance of 
100 pc. The dash - dot curve, dashed curve, and solid curve show stages of evolution at t = 5 x 10 6 yr, 5 x 10 7 yr, and 7 x 10 7 yr respectively. 



Fig. 10.— The spectral energy distribution for a CB disk for the model sequence with a = 0.005 and S = 0.001 for an assumed distance of 
100 pc. The dash - dot curve, dashed curve, and solid curve show stages of evolution at t = 8 x 10 6 yr, 8 x 10 7 yr, and 9 x 10 7 yr respectively. 
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Fig. 11.— The spectral energy distribution for a CB disk for the model sequence with a = 0.005 and 5 = 0.005 for an assumed distance of 
100 pc. The dash - dot curve, dashed curve, and solid curve show stages of evolution at 10 6 yr, 5 X 10 6 yr, and 9 x 10 6 yr respectively. 
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Fig. 12.— The spectral energy distribution for a cataclysmic variable system with CB disk, for a = 0.005 and S = 0.005 at an evolutionary 
time of 5 x 10 6 yr. At this time, M 2 , is 1.3 x 10 18 g s —1 . The distance is assumed equal to 100 pc. The dash - dot curve shows the contribution 
from a low mass secondary star of mass 0.36 Mq, the dashed curve corresponds to the sum of the contribution from a steady state accretion 
disk surrounding the white dwarf component of mass 1.14 Mq and the low mass secondary, and the solid curve includes the contribution from 
the CB disk as well. Note that the CB disk dominates the spectral energy distribution for A > 3/im. 





